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ABSTRACT 



In this paper we analyze one-matrix models by means of the associated discrete 
linear systems. We see that the consistency conditions of the discrete linear system 
lead to the Virasoro constraints. The linear system is endowed with gauge 



invariances. We show that invariance under time-independent gauge trans- 
formations entails the integrability of the model, while the double scaling limit is 
connected with a time-dependent gauge transformation. We derive the continuum 
version of the discrete linear system, we prove that the partition function is actu- 
ally the r-function of the KdV hierarchy and that the linear system completely 
determines the Virasoro constraints. 
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1. Introduction 



In the study of two dimensional quantum gravity, one of the main approaches 
is provided by the random matrix models. In this context the most interesting 
and calculable model is the one Hermitian matrix model, which has been studied 
intensively [1] (for reviews, see [2]). At the discrete level, it is an integrable system 
[3], and its partition function satisfies certain constraints, which belong to the Borel 
subalgebra of the c=l non-twisted Virasoro algebra[4]. After taking the double 
scaling limit the system is expected to become a KdV system restricted by twisted 
Virasoro constraints [5,6]. However, there is still a gap between the discrete level 
and its continuum version: first, the meaning of the Virasoro constraints is not 
clear; secondly, the identification of the partition function with the r-function as 
well as the appearance of the twisting has not been justified yet. 

In this paper we will try to shed light upon these points by analyzing an 
auxiliary tool in one-matrix models: the relevant associated discrete linear systems. 
We will show in fact, in section 2, that a one-matrix model is equivalent to a certain 
discrete linear system (DLS) in which the string equation appears as a compatibility 
condition. In section 3 we will derive the discrete Virasoro conditions from the 
consistency conditions of the DLS. In section 4 we examine the "time-independent" 
gauge symmetry of the DLS and show that it is equivalent to the integrability of 
the system. In section 5 another kind of gauge symmetry is considered which 
includes in particular rescalings of the couplings, and leads naturally to the double 
scaling limit. In the last two sections we examine the continuum limit of the DLS 
and study its properties. In section 6 we determine the continuum analog of the 
DLS and find its consistency conditions. Finally, using the latter exactly as in the 
discrete case, in section 7 we determine the continuum Virasoro constraints. Our 
result is that the partition function is actually described by a KdV-r function, and 
the constraints are the twisted ones. This result had only been conjectured up to 
now. 
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2. One-Matrix Model 

In this section we review the main results concerning the one-matrix model 
and introduce the associated discrete linear system (DLS). 

The partition function of the one- matrix model is defined by [7] 

N N 



i=l i=l 



Z N (t) = J dMe~ TrV ^ = ynrfA l A 2 (A)exp(-^\/(A. t )) (2.1) 
with the potential 



oo 

r 



V(M) = ^t r M 



r=l 



A powerful tool to study this model is the use of orthogonal polynomials. They 
are normalized as follows 

P n (A) = \ n + . . . (2.2) 
and satisfy the orthogonal and recursion relations 

J d\e~ v ^P n (\)P m (\) = 5 n!m h n (t) (2.3) 

AP n (A) = P n+ i(A) + S n P n (\) + RnPn-i(X) (2.4) 
It is convenient to introduce another set of the polynomials 

Vhn 

Then eqs. (2.3) and (2.4) become 

J d\£ n (t, \)£ m (t, =Sn,m (2.5a) 



2 



A£ n = \/Rn+lin+l + S n ^ n + V^Cn-l (2.5b) 

where the function R n (t) are defined as the ratio 

D — ^ n \ 1 

itn = 7 , n > 1. 

One can show that 

d 

S n = - 1 r-lnhn(t), n>0 (2.6a) 

or more generally, 

■^\nh n (t) = -(Q r ) nn n>0 (2.6b) 

where Q is the Jacobi matrix 

Qnm = J d\£ m (t, A)A£ n (£, A) 

= \/ Rn+\$n,m— 1 "I" S n 8 n ^ m + "\J Rn$n,m+\ 

In the same way, we introduce the matrix P 

rfAU(A)— £n(A) (2.8) 



(2.7) 



Using the orthogonal polynomials, we can perform the integrations in (2.1) and 
obtain 

Z N (t)=N\hohih2...h N -i (2.9a) 

<9 iV ~ 1 

\nZ N (t) = -J2Qnn = -TrQ r , r>\ (2.9b) 

Next we introduce the discrete linear system alluded to in the introduction. Let 
us denote by £ the column vector with components £o, £l, £2, • • ■ , use the recursion 
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relations and differentiate the orthogonality relations with respect to t r : we arrive 
at the following discrete linear system (DLS) of equations 



(2.10a) 



(2.10b) 
(2.11a) 



P = Y J kt k Q k a - 1 



(2.11b) 



k=2 



where the dependence on t and A has been understood. Here and throughout the 
paper we adopt the notation 



' \{Q r )nm, m<n 



(Qa)nm — \ 0| 



m=n 



(2.12) 



k -\{Q r )nm, m> n 



The product in the above equations is of course the matrix product. 

The consistency conditions for this linear system give rise to the discrete KdV 
hierarchy [8] 



(2.13) 



and to the so-called string equation 



[Q,P] = i 



(2.14) 



All the integrability and criticality properties of the matrix model are encoded 
in the DLS. Showing this is the aim of our paper. 
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It is interesting to answer the question: to what extent is the correspondence 
between discrete linear systems and one-matrix models one to one? There cer- 
tainly are linear systems that do not correspond to matrix models, however if we 
impose the matrix Q to have the Jacobi form (2.7), the correspondence is one to 
one. Indeed let us start from the infinite column vector £ with orthonormalized 
components £o,£i,£2, ••• as m eq.(2.5a), and write the system 

Qi = K (2.10a') 



= Qr ^ (2 - 10b,) 

Then we can reconstruct the partition function from 



If we define now 



^ r \nZ N (t) = -TrQ r , r>l 



we have the consistency condition 



9 ^ = [Qa,P\ (2-lla) 



r 



This admit the only solution 

oo 

P = J2 kt kQa~ 1 (2.Hb) 

using simply eqs.(2.10a') and (2.10b') beside the orthonormality conditions. 

To end this section let us pause a bit to make a comment on the DLS. As we 
will see, it is justified to consider eq. (2.10a) as a discrete version of the Schrodinger 
equation where A plays the role of the spectral parameter; while eq.(2.10b) shows 
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the KdV-type of (isospectral) deformations of the discrete Schrodinger equation 
we will be discussing later on. It is well known that, typically, the Liouville theory 
is characterized by a Schrodinger equation (actually by two of them, one for each 
chirality). Even though we will be using throughout the paper only the linear 
system (2.10-11), it is interesting to notice that we can push this analogy further 
by introducing a discrete version of the Drinfeld-Sokolov linear system. This is 
done as follows. 

Introduce the two matrices Q\ and Qi = (QiY, t meaning transposition, by 



(Ql)nm — V X n 5 n ^ m + <| <5n,m+l (2.15) 

y A n _i 

and the block matrix and vector 

•-£;)■ -a) 

Then the discrete Drinfeld-Sokolov linear system is 

QE = (2.17) 

This implies in particular = A£i, and we recover eq.(2.10a), provided Q = 

QiQ\- The latter condition allows us to uniquely determine the X n 's introduced 
above in terms of the 5 n 's and i? n 's. One finds 

X = S , X n = ^-, n>\ (2.18) 

where, introducing for the sake of homogeneity the notation E4 = Ru-i, 
Y = S 

Y\ = SiSq — Rio 



Yn - ^2 ( -1 ) fc Sn- • -Sik+lRikik-lSik-2- • ■ 

0<fe<§ ii<i 2 -K...<ifc-l 

' ' ' Si 2 +iRi 2 i 2 -iSi 2 -2 ■ ■ ■ Sh+lRiiii-lSii-2 ■ ■ ■ Sq 



(2.19) 
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We can go further. Let us introduce the sfo generators 

(\ o \ /oi\ /oo 

*=(<> -ij- ^=(0 oj' *-=(i o 

while 1 will denote the identity 2x2 matrix. Next define 

Q1 + Q2 
2 

2 

5+ = E+ + A£_, „4 = n-£+ 
Then the discrete Drinfeld-Sokolov linear system can be written 

(A + ^)S = (2.17') 

We notice that £ + is the sum of the step operators corresponding to the simple 
roots of the affine s/2 algebra, provided we identify the spectral parameter with 
the loop parameter. In this sense A can be thought of as the discrete analogue of 
an SI2 loop algebra connection. 

3. VlRASORO CONSTRAINTS FROM DLS 

An important piece of information for the matrix model is contained in the 
so-called Virasoro constraints' 4 ' 5 ' 6 ^ 

L n Z N {t) = 0, n>-\ (3.1a) 

where 

Ln = J2 kt kjr + 2iV— + 7^7 + N 2 5 n0 , (3.1b) 

[L n , L m ] = {n- m)L n+m . (3.1c) 
They completely determine the possible perturbations. 
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In this section we show that the Virasoro constraints result from the consistency 
conditions of the linear system (2.10-11), eqs.(2.13-14). 

To this end we rewrite the string equation (2.14) in the following form 

fc=2 K 1 

where we have used the KdV equations (2.13). Eq. (3.2) implies that 

00 ft 

eS n = -l, n>0; £ = J2kt k - (3.3) 

which, by (2.6a), can be re-expressed as 

d 

£—\nh n = l, Vn>0 
at 1 

since the operator £ commutes with After integrating over t\, and using the 
formula (2.6b), we get 

00 

J2 kt k(Q k ~ 1 )nn + a = 0, Vn>0. (3.4) 
k=i 

At first sight the integration constant a seems to depend on t2, t%, but using the 
discrete KdV hierarchy and the string equation one can prove that a is actually a 
constant. After summation over n, from (2.9b) we obtain 

{£ - N(h - a))Z N (t) = 

We see that we can absorb a by a redefinition of t\. So finally we get 

00 d 

£^tt Nt 1 )Z N (t) = (3.5) 

which is nothing but the L_i constraint. 
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We remark that choosing even potentials would imply S n = 0; therefore eq.(3.3) 
would be meaningless and would forbid us to recover the L_i Virasoro condition. 
We will see later on that in the continuum limit this obstruction is removed. 

In order to derive the other Virasoro constraints, we introduce the quantities 

EE y/Rn^ . . . R n - r P n , n - r , T > 0. (3.6) 

Due to the string equation (2.14), one finds for the above symbols the recursion 
relations 

E}(H-1) _ R (H-1) _ (q _ c \ R ( r ) , p R ( r— 1) 

J -'n+l J-'n — \°n—r ^>n)- u n ~r rt n — r u n 



The first few ones are as follows 



= 0, B^ = n 



B { n ] = S + Si + . . . + S n _ 2 - (n - l)S n -! 

n— 3 n— 4 n— 4 

-Bn 3 ^ = ^ <Sf + SiSi+l - Si(S n -2 + S n -i) 
i=0 i=0 i=0 

n— 3 



— S n sS n -i + 2 ^ Ri — (n — 2)i? n _2 



(3.7) 



On the other hand, from the KdV equations, it is easy to see that 

Q Q N-l 

w ^=arX Si = - RN - 1 - (38) 

i=0 

Furthermore, since Q r is a symmetric matrix, from the eqs.(2.11), (3.4) with a — 0, 



one finds that 



( Pnh n > I 



Jk-i) 



= < 



0, n=l . 



E kt ^nl 

V -P nU n < I 
Therefore, for any positive integer r, one obtains 



(3.9) 



J2 kt kQnn +r) = J2 E (kt k Q( k -V) nl (Q r +l) 



oo n+r+1 



in 



k=l 



k=l l=n—r—l 



(3.10) 



B (r+1) , B (r+1) , 
n n+r+l + n n + 



Then, after summation over n, the above equations give the L r -constraint. For 
instance, for the first three cases, it is easy to check that 



ktkQnn ~ ktkQ n: n-iQn-l,n + ^kQ^m^Qnn + ktkQ^n+iQn+l, 
k=l k=\ k=l k=\ 



B n 1] + B < n 1 l 1 =2n+1 



E ^kQnn 1 = B%{S n + 5 n+ i) + B n l \s n + 5 n _i) + B n 2) + 5 (2) 



n+2 



= 2[S + S 1 + ... + S n - 1 + (n + l)S ri 



E ktkQnt 2 = B% + + S n 2) (5 n + S n ^ + S n _ 2 ) + S^ 2 (S- n + S n+1 + S n+2 ) 



(2) 



k=l 



Bn\Rn + -Rn-1 + Rn~2 + 'S'n-l + S n Sn-l + S%) 



+ B^^Rn + i? n -i + i?n+l + 'S'n+l + S n S n+ i + S^) 



which after summing over n, and making use of (2.9b), (3.7) and (3.8), become the 
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Virasoro constraints 



L n Z N {t) = 0, 



n = 0,1,2 



(3.11a) 



L n — y j kt)~ 
fc=i 



OO 



5 




(3.11b) 



The Virasoro algebraic structure (3.1c) ensures that the higher order constraints 
are also true. 

4. Gauge Symmetry and Integrability 

As we anticipated in the introduction the discrete linear system (2.10) is char- 
acterized by gauge symmetries which, on the one hand, ensure integrability and, 
on the other hand, allow us to envisage the double scaling limit as a singular gauge 
transformation. This section and the following one are devoted to studying the 
properties of these gauge tranformations. 

Let us consider the following transformation (at fixed t^s) 



where G is a unitary matrix. If the transformed Jacobi matrix Q has the same 
structure as Q, i.e. only the diagonal line and the first two off-diagonal lines are 
non-zero, and, moreover, if 




(4.1) 




(4.2) 



then, we say that our linear system is gauge invariant. 
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Let us examine these transformations more closely by considering the infinites- 
imal transformation 



G= l + eg. 

Then, the invariance requires the matrix g to satisfy the equations 

Q = Q + e[Q,g] (4.3a) 



^g = [Q r a ,g]-[Q r ,g]a- (4.3b) 



A non-trivial solution is 



<7 = Z>*<& ( 44 ) 
k 

where bk's are time-independent constants. By abuse of language we will call this 
a "time-independent gauge transformation". 

Let us consider the case when only is nonzero. Then 



5Q = q_q = £bk[Qj Qkj = -eb k —Q. (4.5) 

otk 

This corresponds to the transformation — > — ebk, which can be rephrased by 
saying that the tuning of the time parameters is realized by means of the gauge 
transformation (4.4). 

This transformation has remarkable properties. On the one hand it can be 
considered as the discrete version of the conformal transformations, on the other 
hand it leads to the integrability of the linear system (and consequently to that of 
the one-matrix model). 
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Let us consider in detail the latter claim. We can think of 5Q given by eq.(4.5) 
as originating from a Poisson bracket in the following sense: 

5Q = e{A r ,Q} (4.6) 

That is 

{A r ,Q} = [Q r a ,Q]. (4.7) 

where A r represents a Hamiltonian to be determined. For the Hamiltonians we 
make the ansatz 

H r = -J2Qnn r = l,2,... (4.8) 

and corresponding to the choice A r = H r , H r -\, iJ r -2, • we obtain different Pois- 
son brackets. More explicitly we write 

{H r - k+1 ,Q} k = [Q r ,Q], (4.9) 

In the following we will study in detail only the cases k = 1,2,3. Comparing the 
LHS with the RHS which can be explicitly calculated, we can obtain the Poisson 
brackets for R{ and Si. Of course we have no a priori guarantee that the brackets 
so obtained satisfy the Jacobi identity. This has to be checked a posteriori. 

While explicitly working out the Poisson brackets one realizes that there are 
two distint regimes according to whether Si — or ^ 0. 

i) First regime, i.e. Si = 0. We find two meaningful Poisson brackets: 

{Ri, Rj}i = R i Rj(5j j i + i - Sij + i) (4.10) 

and 

{Ri, Rj}s =R{Rj(Ri + Rj)(5j ; i + i — <5i,j+i) 
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+ RjRj-iRj-28j,i+2 — RiRi-lRi~2$i,j+2 



(4.11) 



while 

{Ri,Rj} 2 = 

ii) Second regime, i.e. Si ^ 0. We have two Poisson brackets: 

{Ri, Rj}i =RiRj(5 jji+1 - 8ij+i) (4.12a) 

{Ri, Sjh =R i S j (5 jyi+1 - 8 id ) (4.126) 

{Si, Sj}i =Ri5j j i + i - RjSij+i. (4, 12c) 

and 

{Ri, Rj} 2 =2R l R j (S i 5 iJ - 1 - SjSij+x) (4.13a) 
Sj}2 =RiRj(b~ij-i + 5ij) — RiRj+i(Sij+\ + Sij+2) 

+ RiS](5ij-5i J+1 ) (4.136) 

{Si, Sj} 2 =(Si + Sj^RjSij-x - RiSij+i) (4.13c) 

For k = 3 eq.(4.9) does not define a consistent bracket. 

Let us conclude this section with a few remarks. Due to eq.(4.7) one would 
expect all the Hamiltonians to commute 

{H n ,H m } = 0. (4.14) 

for any meaningful Poisson structure. However due to the subtleties connected 
with traces of infinite matrices, one should verify this property starting from the 
Poisson brackets (4.9) and (4.10). We have done it for the first few cases. 
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Moreover, from the definition of the Poisson brackets, we know that 

{H r+2 , Q}\ = {H r , Q}s 

in the first regime and 

{H r+1 ,Q} 1 = {H r ,Q} 2 

in the second, namely the two Poisson structures are compatible with each other. 
That is to say, in both regimes the linear system has infinitely many conserved 
quantities and possesses a bi-Hamiltonian structure. 

Finally the Poisson brackets (4.10) and (4.11) are the same we come across 
in the lattice version of the Liouville model [9]. In that case Ri plays the role 
of a lattice deformation of the classical Virasoro generators. However one should 
remember that in the first regime we cannot impose the string equation (see remark 
after eq.(3.5)). So we find a remarkable similarity of structures in the Liouville 
model on the lattice and in the one-matrix model unconstrained by the string 
equation. Strictly speaking, the Poisson structures of the one-matrix model are 
those of the second regime. It is an interesting open problem what field theory on 
the lattice they correspond to. 

5. Reparametrization and Time-Dependent Gauge Transformations 



The DLS (2.10) is form invariant under reparametrization of the couplings. 
That is 

Q(tK(i) = -JU(t) = (Q r a )(tm 

Ot r 

where t = (t\,t2, •••) and tk is a smooth functions of the t^s. 

This invariance is a formal one. However, by combining gauge and reparametriza-| 
tion transformations, we can obtain significant symmetries of the system. Let us 
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consider the transformations 



tk — ~tk + £ (k ~ n )tk-n, Vk > n, n > — 1 
t n = t n - 2Ne, n > 1 

[Q(t) = G- 1 (t)Q(t)G(t) 



g = P, n = -l 

g = QP, n = 



(5.1) 



In this kind of setup it is possible to find G = 1 + eg so that 

Q(t) = Q(t) (5.2) 

and the linear system becomes 

' Q(t)i(t) = xi(t) 

i-J(t) = Q r a(t)i(t) 

We refer to these as time-dependent gauge transformations. 
Let us consider two examples 



(5.3) 



(5.4) 



In these two cases eq.(5.2) is satisfied (the linear system is invariant but remark that 
£ is not invariant for n — 0) and this fact leads straightforwardly to the L_i and 
Lq Virasoro constraints. We could as well obtain the other Virasoro constraints, 
but in these cases the matrix g has a complicated form and will not be written 
down here. 
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6. The Double Scaling Limit 



The purpose of this section is to recover a continuum version of the DLS (2.10) 
in the double scaling limit. In general we will exploit the idea that the double 
scaling limit is mimicked by a finite version of the transformation (5.1) above 
when n = 0: 



where 7 is a finite constant. We recall that under this transformation Q remains 
invariant. In other words, the double scaling limit is connected with a singular 
case of a symmetry operation on our DLS. 

Let us consider a k-th order critical point and define, as usual, the continuum 
variables 




Vr 



(6.1) 



n 



R(x) = R, 




Moreover we set 



e = 




*o = (l 





(6.2) 



The double scaling limit corresponds to 



(3 -> 00, 



N 



00, 



to fixed 



For large n ~ iV one has 



x = 1 - e 2k t , 



R(x) = 1 + e 2 u(t ) 



(6.3) 



where u(to) is the specific heat. 
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The latter ansatz requires a comment. Let us consider the string equation 
(2.14), suitably rescaled 



. oo 
[Q,P] = 0, P = J2 r trQ r a~\ t r = Pt r (6.4) 



r=2 



where t r are renormalized coupling constants (see below). This equation estab- 
lishes strong restrictions between the limiting expressions of Q and P. With the 
simplifying assumption t2r+i = Vr, the above ansatz is correct, as is well known; 
if we switch on the odd interactions the analysis is more complicated, the above 
ansatz is still correct but we have not been able to exclude other solutions. In the 
following we will stick to the case of even potentials and to (6.3). 

Let us now write down a few expansions which will be useful in the following. 

It is easy to see that 

R nT i = R(x =F — ) = 1 + e 2 u(t ± e) 



Then, using Taylor expansion, we have the following formulas 

2 o ' b 2 a a z 
Rn+b = l + ew — be u + — e u — — e u + . . . (6.5a) 



1 1 

Y[R n+b+i = l + {q + l)e 2 u ~-{q + 1)(26 + q)e 3 u 



i=0 



+ t *(±^ + (l±plA (8.5b) 



i=0 
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A rs _ i , 9 + 1 2„ (g + l)(26 + g) 3 ./ 

XX V u n+b+i - 1 H ^— e u ^ e u 



i=0 



+ A±(t^l u ' + (i±m^iiA (6 . 5c) 

_ e , f f (±±3l u -> + (» + »)(« + Dfa-D ^ + . . . 
V ' 12 8 / 

Using these formulas one can derive 

Q = 2 + e 2 (d 2 + u) + 0(e 3 ), (6.6a) 

-Q a = ed+ le 3 (d 3 + 3ud + ^u) + \e\2ud + u) 
6 2 8 

1 * ' 1 2 -}3 , 'a2 



+ -e° —9° + -m<9 j + + « 9 - m 2 <9 - mm + -« + • • • , (6.66) 
8^15 3 3 J 

-Ql = 2ed + i 3 (d 3 + Ivd + \u) + \e\2ud + u) 

c/4r 4 o ' n 3 " 5 "A 
+ e 5 — d 5 + -m<9 3 + 2m <9 2 + -m ,9 + — m +..., (6.6c) 

V lo o & LA J 

-Ql = 6ed + e 3 (5<9 3 + 9ud + \u) + V(2m'<9 + u) + ^ (^<9 5 

,„ _o 45 ,-9 37 «. 9 9a 9 / 11 »A / n n 7\ 

+ 15m<9 3 H u <9 2 H it d + -m 2 <9 + -mm H u )+..., (6.6d) 

2 2 2 2 2 J 

-Ql = 12ed + 16e 3 (<9 3 + ^ud + ^u') + 6e 4 (2ud + u") + 2e 5 (^<9 5 

+ 16m<9 3 + 24m d 2 + 19m d + 6m 2 <9 + 6mm + —u )+..., (6.6e) 
-QS = 30e<9 + e 3 (45<9 3 + 75m<9 + + 7be A {2ud + u") + -e 5 (29<9 5 
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O ' O '' O I III 

+ 90m<9 3 + 135m d 2 + 111m d + 45m 2 <9 + 45mm + 33m ) + . . . (6.6/) 



etc. 



Let us see now some consequences of the above expansions. First of all let 
us notice that in the continuum limit the reduction to even potentials does not 
contradict the string equation as in the discrete case. We should remember that 
the contradiction is exposed in eq.(3.3). From the above expansions it is not 
difficult to see that in the continuum limit it does not make sense to single out an 
equation like (3.3), while the LHS of the string equation is replaced by a differential 
operator even if t2r+i = Vr. So in the continuum limit on can safely choose an 
even potential. 

Next we consider the continuum limit of (2.10a). From eq.(6.6a) we see that in a 
neighbourhood of the critical point A ~ 2. Therefore we introduce the renormalized 
quantities 

A = e- 2 (A-2), Q = d 2 + u (6.7a) 

-if = P£, P = e 2 P = e l ~ 2k P. (6.7b) 

d\ 

Then, the discrete Schrodinger equation goes over to its continuum version 

(d 2 + u)i(t )=\i(t ) (6.8) 
Similarly the string equation becomes 

[Q,P] = 1. (6.9) 



We are now in a condition to explicitly determine critical points. From eq.(2.11b)| 
and (6.7b) (recall that we are working with the simplifying assumption of even po- 
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tential) we have 



P = e 2 (2t 2 Q a + ±UQl + Qt 6 Ql + . . .) 
= e l ~ 2k {2i 2 Qa + 4hQl + Qt 6 Ql + ■■■) 



(6.10) 



We remarked above that the string equation (6.9) puts severe restrictions not only 
in the discrete case but also in the double scaling limit. From eq.(6.6) and (6.7), 
we see in particular that all the operators Q^s are vanishing in the limit e — + 0, so 
that if one wants the string equation to be satisfied, one must let a certain subset 
of bare coupling constants in P go to infinity (DSL). The practical recipe is to look 
for combinations of i 2 rS such that all the singular terms in the second expression 
of (6.10) vanish. Let us see a few significant examples. 

i) k=2. In this case only t 2 and £4 are nonzero, and (3 = e~ 5 . Then, from 
(6.6b,d) we see that only if we set 



are we able to eliminate the e d term in P, and we get the known operator 



As is well-known at the critical point t\ = the above string equation is the 




P = h 2 (d 2 + u)i + 0(e) 



The string equation becomes 



-t 2 n' 2 [u] = i 



where we have introduced the Gelfand-Dickii polynomials 



TZ k [u] = [(d 2 + u) k + Kd 2 + u]. 



(6.11) 
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Painleve equation of first kind 



1 " 2 

^o=g^ +u . 

ii) k=3. Only £2,^4 and t$ are non-vanishing, (3 = e~ 7 . According to the above 
recipe we kill all the negative powers of e in P if we put 

+ 105-71 c -77 

It is straightforward to show that 

P= 7 -Ud 2 + u)\ + 0{t) 
and the string equation becomes 

- 7 -m,[u] = i. 

The third critical point is at t§ = — ^| and the differential equation is 

. O 1 /2 // 1 (A\ 

to = -(u - -u -uu + —u y i). 



One can proceed in this way and determine higher order critical points. As 
is well-known, on a very general ground the form of the operator P must be as 
follows 

00 Ii 

p = J2( n+ 2 )in ®+ + ~ 2 + ° {e) - (6 - 12) 

n=l 

We conjecture that this form is induced by the following coupling redefinitions: 

oo 1 oo 

t 2r = + ^tnCna^e-^+V = ^1% (6.13) 

n=r n—r 
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where 



>) _ ( iy+i 
r ~ [ ' (n-r)\(2r)\ 

We checked eq.(6.13) for the first few cases and found 

Ci = l, c 2 = ~ c 3 = | c 4 = ~... 

In general one has 

c _ ( n n+i (2n-l)!! 

These factors are determined in such a way as to reproduce the standard KdV 
hierarchy. 

It is worth noticing that 1) in eq.(6.12) and (6.13) we are considering all 
the critical points at a time, and 2) the time transformation (6.13) is made of 
a reparametrization plus a scale transformation of the type (6.1). 

What is left for us to do is to analyze the continuum limit of the KdV hierarchy. 
On the basis of eq.(6.13) one naively has 

JL = _ y {n + I )Cne -(2-+i) a M JL. (6.i4) 

dt n ^ 2 ; n dt 2r 1 ^ 

So, in particular, on the basis of (2.10) and (6.6) we must have 



|e=(|^+(9 2 + 4 + o (e )^ 



dt 3 

etc. These are however naive formulae since {t n , n > 1} is not a complete set 
of parameters after we take the continuum limit. To correct this we have to allow 
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also for a ^-dependent term in the RHS of (6.14). This additional term takes 
care exactly of the divergent terms (in the e — > limit) in the RHS of the above 
equations. 

Finally the evolution equations become 

J-f=(d 2 + u)^f, n >0 (6.15) 
dt n 

which result in the standard KdV flow 

J-m= [{d 2 + u) n + +1 \d 2 + u], n>0 (6.16) 

In eq.(6.15) £ is the limit of £ possibly multiplied by an e-dependent factor 
which may be necessary in order to obtain a finite result. 

7. The Virasoro Constraints in the Continuum Limit 

In the previous section starting from the DLS we have obtained, near criticality, 
a continuous linear system 

(9 2 + u)|=A| (7.1a) 

9 z /o2 , sn+i 
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11 

^ = E( n+ oWr 5 ( 7 - 4 ) 



n=l 

whose consistency conditions are 



[Q,P] = 1 (7.2a) 

J-„=[(d 2 + „)^,S 2 + „] (7.26) 

0t n 
0t n 



P=[(d 2 + u) n + + \P}. (7.2c) 



We want now to recover the Virasoro constraints in this continuous system. 
The strategy is the same as for the discrete case. We use essentially the string 
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equation (7.2a). First of all, as is well known, in the continuum limit the partition 
function behaves like 



N-l 

In Z = (N - k) In R k + constant terms 
k=l 
to 

/~J ~ ~J ~J 

dt (to — t )u(t ) + (9(e) + regular terms 



d 2 \nZ = u(t ) (7.3) 

where, in taking the continuum limit, we passed to the normalized partition func- 
tion 

Z = Z(t)/Z(T) (7.4) 

the parameter T being a reference point connected with the extremum of integra- 
tion to — 0. 

Now, using (7.1c), eq.(7.2a) can be written 

-E( fc+ 9)**^*M = 1- (7-5) 



T 
k=i 



Integrating once with respect to to, we obtain 



OO | 

X> + 2^(K k [u] - K k [0]) + t = (7.6) 



fc=i 



i /t-« r^i • /t-« r i i,r r>^7 'onedirectly -consider sthecontinuumlimitofthediscreteVirasoroconstraintsAtisni 

where 7tfc|0J is 7cfc[u] computed at to — 



hisisactuallythecase. . 
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In order to simplify the next formulas let us introduce the recursion operator 



4> = -d 2 + u + -u'd- 1 

4 2 

and define the recursion relation for the Gelfand-Dickii polynomials 

K n+1 = &l n = 4> n du. (7.7) 
Remembering that, on the basis of our conventions, we have 

d- 1 n' k [u] = n k [u]-n k [o] 

eq.(7.6) can be rewritten 

00 -. 

F = t + J2( k + Ad^^du = 0. (7.8a) 
k=l 

On the same basis we can write 

d~ l j) n+l dF = 0, n > -1 (7.8b) 

To obtain these equations we have used only the string equation. We can as 
well envisage eqs.(7.8) as a consequence of a symmetry of the system, precisely as 
a consequence of the fact that u(to) and the KdV hierarchy are invariant under 
the transformations 

h — > t k = i k + t{k - n + \)tk- n 

(7.9) 

u(t") —<£ n+1 •! + «(**) 
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The generators associated with (7.9) are 



00 1 d 1 

£ - 1 = B*+5fc s ^ + ^. 



00 19 1 

k — 



Here, for later purposes, we have introduced a constant p (for example, by rescaling 
all the f's). 

Even though we will not use it in the following, it is worth mentioning that 
there is a larger symmetry of the system: the latter is also invariant under the 
transformations 

+ £ (7.11a) 

whose generators are given by 

\4 = J-, k>0. (7.11b) 

dt k 

The generators V k s and L n 's characterize the master symmetry of the KdV hier- 
archy [9,11]. The corresponding algebra is 

\V k ,Vi] = 0, k,l>0, (7.12a) 

[V k ,L n ] = (k + ^)V k+n , k>0,k + n>0, (7.126) 

[Vb,L_i] = -U (7.12c) 
2p 

[L n , L m \ = (n - m)L n+m , n,m> -1. (7.12d) 
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But let us return to the derivation of the Virasoro constraint. Representing 
now (7.8) in terms of the partition function, we get for n = —1,0 



(00 r\ \ 



or, in general, 

3^/^=0, n>-l. (7.13) 

where, by definition, 

k = Lq - — , l n = L n , n^O. 
16 

The most general solution of (4.28) has the form 

/ n V / I = 6 n V / I (7.14) 

where the 6 n 's are tcrindependent but arbitrary functions of the other parameters. 
In order to determine them we remark that they must be compatible with the 
algebra (7.12). In particular they must satisfy 

[In, bm] - [l m , b n \ = (n - m)b n+m + ^n5 n+ m,Q- (7.15) 

Moreover we remark that (7.12) is a graded algebra provided we define the degree 
as follows 

d 

deg[t fc ] = —k, de S^] = Z ' de sH = 1 
Therefore, deg[L n ] = n and, from (7.14), 

deg[6 n ] = n, n > -1 (7.16) 
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The general form of b n will be a sum of monomials of the following type 

M n (p,qi,..., q a ) = const (? t q n \...t q n a a 

where p is a real number and qi,q2, ■■■ are nonnegative real numbers (we exclude 
negative exponents as it is natural to require a smooth limit of b n as any one of 
the couplings vanishes). Next we remember that the parameter p appeared on 
the scene because of a rescaling t n — > ^/pt n . Therefore if we perform the opposite 
rescaling the dependence of b n on p must disappear. This implies the condition 
V = \{qi + ••• + Qa)- So the degree of the above monomial would be 

a 1 

deg[M n (p, 91, Qa)] = J^(- - ni)qi (7.17) 



Comparing eq.(7.16) with (7.17), we see all the 6 n 's are zero except perhaps for bo 
which must be a constant, and b-\, which could depend on t\ and p. We now use 
the consistency conditon (7.15). For n = 1 and m — — 1 it tells us that 

5~ <%_i OJ 1 
-t 2 -^- = 2b + - 
2 dh 8 

This allows us to conclude that 

b0 = -1. 

16 

and b-\ does not depend on t\. Next we use again (7.15) for n = and m — — 1 
and conclude that 

6_i = 

Collecting the above results we obtain the Virasoro constraints 

L n v / I = 0, n>-l (7.18) 

with L n given by (7.10). Another proof of the same result is given in Appendix. 
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Finally we recall that due to (7.2b) and (7.3), \/~Z is a r-function of the KdV 
hierarchy. 



Appendix 

This Appendix is devoted to another proof of eq.(7.18). We start from eq.(7.15) 
for m = — 1 and m = and analyze the to dependence. Since all the b n are Iq- 
independent we get, in the first case (m = —1), 

—J- = 0, n > 1 

ot n 

and, in the second case (m = 0), 

^ = 0, n>l 

Ot n 

In conclusion both bo and 6_i are constant. Applying now again eq.(7.15) for n = 
and m = — 1, we obtain 

= 

i.e. eq.(7.18) is true for n = —1. Since we have shown above that is a 
r-function of the KdV hierarchy, we can now apply a theorem of ref.[ll] which 
asserts that if L-\t = 0, then L n r = 0, Vn > —1. Consequently eq.(7.18) is 
proven. 
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